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Abstract. We prove that in the Hartle-Hawking approach to quantum cosmology the 
existence of an inflationary phase is a general property of minisuperspace models given 
by a closed Friedmann-Robertson- Walker universe containing a massless scalar fleld with 
a Xcf)^ self-interaction. The evolution in time of the cosmic scale factor and of the scalar 
fleld in the very early universe is derived, together with the conditions to be satisfled in 
order to solve the horizon and flatness problems. 
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1. Introduction 

The quantum state of the universe [1, 2] is described by a functional ip of the 3-geometry 
and of the matter field configuration on a compact spacelike 3-surface S. Any solution 
of the Wheeler-DeWitt equation and of the momentum constraints is a possible quantum 
state of the universe [1-3] . In the last few years there has been great interest in the Hartle- 
Hawking proposal for the boundary conditions of the universe [1, 2, 4, 5], which seems 
to predict a universe whose large-scale features are in agreement with observations [1, 
6-9]. According to Hartle and Hawking it seems that ■0 is to be taken over all compact 
Euclidean 4-metrics which induce the given 3-metric on their boundary S, and over all 
field configurations which match the given value on S and are regular on the compact 
4-manifolds having S as their only boundary. 

In order to solve the Wheeler-DeWitt equation, most authors have considered the 
minisuperspace approximation [10]. In the by now well known Hawking massive scalar 
field model [1, 2, 11, 12], the universe begins its inflationary expansion from a non-singular 
state. At the end of the inflation, it undergoes a matter- or radiation-dominated expansion, 
reaches a maximum radius and recoUapses to a singularity. However, in general it is not 
clear whether the solutions to the classical field equations corresponding to the Hartle- 
Hawking wavefunction describe a universe which will recoUapse to a singularity. In fact, 
Louko [13] has pointed out that the method of providing initial conditions for purely 
Lorentzian solutions by analytically continuing purely Euclidean solutions in the matter- 
dominated and recoUapsing phases may be erroneous. His calculations seem to suggest 
that the universe will not reach the singularity in the future. 

After this short exposition of some ideas and problems of quantum cosmology, let 
us now turn our attention to the problem we try to solve in this paper. We study a 
minisuperspace model given by a closed Friedmann-Robertson- Walker (hereafter referred 
to as FRW) universe containing a massless scalar field with a Xcj)'^ self-interaction. A 
massless scalar field is considered for two main reasons : 
(a) it is not clear whether a fundamental massive scalar field exists [7] ; 
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(b) one is still unable to say which theory of particle physics gives the best description of 
the very early universe [14]. 

In section 2 we write down the action integral and the field equations of our minisu- 
perspace model, and we discuss in detail how the Hartle-Hawking proposal enables one to 
impose initial conditions for the solutions to these equations. We also endeavour to derive 
the oscillatory region for the wave function in view of its behaviour in the limit of small 
3-geometry and by comparison with the massive scalar field model. In addition we discuss 
how the spacelike or timelike nature of the surfaces of constant potential of the Wheeler-De 
Witt equation may be useful in deriving the behaviour of the wavefunction. 

In section 3 at first we focus our attention on the region in which the wavefunction 
oscillates. The phase S of the WKB approximation is the one obtained by analytical 
continuation of the Euclidean action for compact metrics and regular matter fields. We 
are thus able to derive the evolution in time of the cosmic scale factor and of the matter 
field using an approximate form assumed by S. The Hartle-Hawking trajectories so derived 
are exact solutions to the Lorentzian field equations. They are distinguished with respect 
to the general set of solutions in that they are inflationary and singularity- free in the past. 

In section 4 we derive the constraints to be satisfied in order to solve the horizon and 
flatness problems; in section 5 we summarize the results obtained and we mention some 
related problems to be studied. 



2. Minisuperspace model and boundary conditions 



Our model is based on the FRW metric which may be locally cast in the form: ds^ = 
a'^{-N'^dt^ + a^{t)dQl), where a = , N is the lapse function [15], a{t) is the cosmic 
scale factor and dOg is the metric on a unit 3-sphere. 

Let (j)'{t) = -^^^ be the scalar field (in these units (j) is dimensionless) having a 

potential of the form V{(j)') = — |— • The action integral containing a surface term [5] 
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becomes, in our case, 



J= -1 / Na' 



q2^2 ^2 ^2 



dt^ Ldt (2.1) 



where A = AV-"+^(7rv^)-^+2. 

Before commencing the calculations, it is worth emphasizing the main property of 
minisuperspace models containing scalar fields. The minisuperspace model given by a 
homogeneous isotropic universe without matter fields but with a positive cosmological 
constant (De Sitter space) did show that the Hartle-Hawking proposal leads to inflation if 
some mechanism is able to give rise to an effective cosmological constant in the early uni- 
verse [1, 11]. But our universe is not expanding exponentially. This is why we require that 
the parameter we call the cosmological constant has finally to vanish as the time passes. 
The simplest model for generating a decaying cosmological constant is the one containing 
a massive scalar field. In fact, when is initially very large and almost constant, the ■m?'(fp' 
term in the action of the massive model acts as an effective cosmological constant. This 
implies that, when =constant=0o7 ah formulae for the massive model can be obtained 
from the De Sitter model by putting H — ^Jrn?(^ — mcpo, where H is the parameter of 
the De Sitter model related to the cosmological constant Ac by means of the relation [1] 
H'^ = If a < H~^, it is known that there are two solutions of the Euclidean field 

equations, and the dominant contribution to the semiclassical approximation of the wave 

function comes from the one having action [1] Ie — ~2jj2 1 — (1 — H'^a^)i , where Ie is 

the action of the smaller part of a 4-sphere of radius ^ bounded by a 3-sphere of radius 
a. In our model, if the term Xcf)^ in (2.1) has to act as an effective cosmological constant 
when (p is large and constant, the parameter H becomes H — -s/XcpQ — ^/X4>q ■ From now 

on, H will always denote VX(f)Q in our equations. 

We can now turn again our attention to (2.1). Putting iV = 1, the Lorentzian field 
equations are 

ad+^ + l + la^^-lxa'r = (2.2) 



4 



Inflationary solutions in quantum cosmology 

+ 3-0+^A(/>^-i =0. (2.3) 

We now propose to derive the initial conditions for the solutions to (2.2) and (2.3) implied 
by the Hartle-Hawking proposal. To this purpose let us recall that the Hartle-Hawking 
ground state is a Euclidean functional integral taken over compact 4-metrics and regular 
matter fields. In the Euclidean regime one thus has [1] 

a(T = 0) = a(r = 0) = 1 (2.4) 

0(r = 0) = 00 0(r = 0) = 0. (2.5) 

In fact, having taken the path integral over compact 4-metrics, a(r) has to vanish at 
some value of r we can choose to be zero [1]. This implies that at large values of one 
has a(r) = -^sin(iyr), and this relation is in agreement with (2.4). In addition, in the 
Euclidean regime (2.3) can be written in a similar form, provided we take the derivatives 
with respect to r and we change the sign in front of ^At/)"^"^. Therefore, multiplying both 
sides of (2.3) by a, and letting a go to zero, one has : a</) ^ 0, ^a0"'~^ 0, in view of 
the regularity of 0. Therefore one has a(3^0) = if and only if = 0, in view of (2.4). 

In quantum cosmology one is mainly interested in the oscillatory regions for the wave- 
function, because such regions correspond to classically allowed regions [1] . One can then 
make the WKB ansatz: ^ = ^ -RefC^ exp(iS'„)] (in what follows we shall consider a single 
component of the wave packet). The phase 5' is chosen to satisfy the Hamilton- Jacobi 
equation for general relativity, and the Hartle-Hawking proposal picks out the Hamilton- 
Jacobi function which is the analytic continuation of the Euclidean action. This implies 
that the Lorentzian Hartle-Hawking trajectories are obtained by analytic continuation of 
the compact Euclidean paths with regular matter fields. The application of this method 
yields, putting t — ^ it [12]: 

ait) = \- cosh{Ht) (2.6) 
H 

m = 4>o (2.7) 
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at very small values of the time. Thus the minimum value of a in the Lorentzian regime 
is equal to the maximum value of a in the Euclidean regime, and the initial conditions for 
a differ greatly from (2.4) because one has a{t = 0) = a{t = 0) = 0. Therefore at large 
initial values of (p, the initial conditions for the solutions to (2.2) and (2.3) are (see also 
(3.3) and the subsequent discussion): 

(f){t = 0) = 00 a{t = 0) = — ^ = ao (2.8) 

^{t = 0) = a(^ = 0) = 0. (2.9) 

Moreover, putting [8] a — log(a), p'^ = — ^ ^^gf ° ^ , the Wheeler-De Witt equation for our 
model becomes 

2 

Now, the action / of a compact solution of the field equations vanishes in the limit 
of small 3-geometries and one often argues that the prefactor A of the semiclassical ap- 
proximation of the wavefunction, ip A cxp(— /) (this is the formula for the solution 
of the Wheeler-DeWitt equation to leading order), may be taken to be a constant, so as 
to normalize to 1 in the limit of small 3-geometries [8, 16]. This would be a boundary 
condition for (2.10) implied by the Hartle-Hawking proposal, and in numerical calculations 
it has been used so as to derive the oscillatory or exponential regions of the wavefunction 
[16]. Thus, from the knowledge of ij; at small 3-geometries, one derives its behaviour at 
larger 3-geometries [16]. 

Strictly speaking, this is not completely true, because the prefactor of the semiclassical 
approximation is given by the path integral of exp(— /2), where I2 is the part of the action 
quadratic in the fluctuations about classical solutions. This calculation is not trivial and 
one can show, for example, that in the case of pure gravity the prefactor diverges in the 
limit of small 3-geometries in the Euclidean regime [17]. Thus, in general, the normalization 
of -0 to 1 may not be sufficiently accurate. However, in our opinion one can still take, to 
a first approximation, the point of view of Hawking and Wu [16]. They point out that 
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a suitable choice (also made in equation (2.10) in this paper) for the factor ordering in 
the Wheeler-DeWitt equation is consistent with a constant behaviour of the prefactor and 
of the wavefunction in the limit of small a, as can also be verified in our case. In fact 
(i^ — vanishes if ■0 is a constant and the potential V{a,(f)) = e^'^Xcf)^ — e^" = 

a^Xcj)"^ — a^, for a given value of 0, vanishes when a — > 0. The numerical integration of the 
Wheeler-DeWitt equation with this condition shows that in the massive model starts 
to oscillate when | | is greater than one and V{a^ 0) is positive. In our case, it is clear 

1 n 

that, even if one puts = m, (j)^ = one cannot cast (2.10) in a form which is formally 
identical to the massive case, because -g^ is very different from However, in the 

case of a theory, for a given value of one has an even larger range of values of a 
for which the potential V(a,0) in the Wheeler-De Witt equation is positive, if A is not 
too small. The sign of this potential plays a crucial role in determining the oscillatory 
behaviour of and this is a quite general result not depending on the particular theory 
for the scalar field [18]. It is particularly useful to discuss this matter in detail. To begin 
with, we point out that the Wheeler-De Witt equation (2.10) is a hyperbolic equation in 
the minisuperspace with metric ds'^ = e^"(— <ia^ -|- (i0^), and we consider the regions in 
which the surfaces of constant V = e^"A0"' — e^'* are spacelike with respect to this metric. 
The idea is that in a local region we can perform a Lorentz transformation [18] to new 
coordinates (5,0) so that the constant-F surfaces are parallel to the axis. This implies 
that (2.10) takes the form — + V{a) ip — and can be solved by separation 

of variables, putting ■0 = '0i(a)'02(0)- One thus finds the differential equation for '01 (5): 

-|- A -|- V{a) tpi{a) — 0, whose solution has oscillatory behaviour if A -|- V{a) > 0, 

where A is a separation constant. It is also possible that the surfaces of constant V are 
timelike. In such a case they can be locally rotated so as to lie parallel to the a axis, which 
in turn implies that V = V{(j)). In such a case one finds the following differential equation 

for '02(0): + A' — V(0) 02(0) = , whose solution has oscillatory behaviour if 



A' — 1/(0) > 0, where A' is another separation constant. Thus V is expected to be negative 
in this case, if the separation constant A' does not assume too large a positive value. 
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To sum up, if the V = constant surfaces are spacelike in the metric of the minisu- 
perspace, ijj is expected to oscillate if V>0 and to behave exponentially if V<0, whereas 
the converse should hold true if the V = constant surfaces are timelike. In our model, the 
nature of the surfaces of constant potential is determined by the quantity E = g°'^V^aVbj 
where g"''' = e^" diag(— 1, 1). The sign of this quantity shows the timelike or spacelike 
nature of the vector normal to these surfaces. Computing S when V = k > 0, one finds 



' " <0 



V=k>0 

2 ■ 



, {k + e*")^ + 2{k + e*")(5/c + e^") + {5k + e^")^ 

0^ 



if (j) or 6°^ are sufficiently large. This implies that the surfaces of constant positive potential 
are spacelike for such values of 0, which in turn implies that the wavefunction has to 
oscillate. Honestly speaking, in the massive scalar field model there are some regions in 
which this method disagrees with numerical results for the wavefunction derived under the 
assumption that ■0 = 1 when a ^ [19]. In view of the fact that both these methods have 
flaws, the question seems to us to be an open one. We can thus show that the wavefunction 
of our model oscillates when V{a, 0) is positive in one of the following two ways: 

(a) by comparison with the results derived in the massive model, pointing out that 
also in our case one can take ^ — > 1 as a 0, and that in our case V{a, (p) has an even 
stronger positivity property if the self-coupling parameter A is not too small; 

(b) by studying the nature of the surfaces V{a, 0)=constant. 

In the case of a theory, we are in agreement with the result derived in a paper by 
Zhuk [20] . However in this paper the normalizability of V' to 1 is not a boundary condition 
assumed from the beginning, but a consequence of the mathematical derivation, in which 
one requires that the solution of the Wheeler-De Witt equation has to vanish when the 
cosmic scale factor oo. 

So far we have emphasized that, at almost constant values of 0, our model can be 
obtained from the massive one replacing the positive quantity m^0^ by means of A0". 
This is possible only when is positive, or when is negative but n is even. In view of 
this result we shall not take into account the case when is negative and n is odd. Of 
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course, the initial value of | | is not completely arbitrary, because the form of S involving 
(f) must be a solution of the Hamilton-Jacobi equation in the oscillatory region for the 
wavef unction. In our case S satisfies the equation 

= e^-- \<j,-e^- (2.11) 



\d(j) J \da 

where the right-hand side of (2.11) is negative. A careful examination of the quantity 
E shows that this is possible only if | (/> | is initially at least greater than a number of 
order one. We have not solved numerically the Wheeler-De Witt equation for our model, 
and this is an interesting problem to be studied for further research. However, it is to be 
remarked that the numerical techniques used so far have some limits. For example, the 
leapfrog algorithm is a valid approximation only when the potential V is everywhere much 
smaller than the inverse of the square of the grid step size [19]. 



3. Lorentzian Hartle-Hawking trajectories and field equations 



In this section we shall perform at first our calculations for positive values of 0, and later 
on we shall also discuss briefly the case when n is even and is negative. 

The Euclidean action for compact metrics and regular matter fields is given in our 



model by Ie 



1 — (1 — e " A0"^) 2 . The geometrical interpretation of Ie is that it is 



the action of the smaller part of a 4-sphere of radius , bounded by a 3-sphere of radius 

a = 6°^; in so doing, we generalize what has been done in the case of the massive model [1, 
11], as already emphasized in section 2. The solution of (2.11) is, in the Hartle-Hawking 
approach, the analytical continuation of Ie- At large values of it is given by 

^ (e2"A0^-l)i. (3.1) 



3A0 
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In fact the insertion of (3.1) into the left-hand side of (2.11) yields 



d(t)J \daj 402 ^ ^ ^ 9A202n+2 



„2 2a 

(e^^A^'^ - 1)2 + e^" - A0"e6" (3.2) 



3A<?i"+2 

and we can easily recognize that the right-hand side of (3.2) reduces to e"^" — A^'^e^'* at 
large initial (j) and small t. The relation (3.1) shows again that (f) cannot be negative when 
n is odd, because in such a case S becomes complex, whereas by its very definition S is 
real. In addition we know that the function S in (2.11) and (3.1) defines a first integral of 
the system: 

_dS _dL 

Pa — — TTT 

oa oa 
_dS _dL 

where L is the Lagrangian defined in (2.1). This system therefore becomes, putting N = 1 

a= (a^Ac/)" - 1)^ (3.3) 

It is to be remarked that this system is in full agreement with the initial conditions 
(2.8) and (2.9) in the Lorentzian regime derived from the Hartle-Hawking proposal. In 
fact, for example, a{t = 0) = 0, when substituted into (3.3), yields oq = — As the 

time passes, the dominant contribution to a{t) is given by the increasing exponential, so 
that a{t) = ^^^2H*^ ■ Under such assumptions, (3.1) takes the approximate form: 



e"^\/A0 



3a ^ 

(3.5) 
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The relation (3.1) enables us to also compute the first correction to (3.5) involving 
negative powers of (f). In fact, from (3.1) one has that 

3 



3A0'^ 



'e2«A0" ( 1 ^ 



6 1 - o , . (3.6) 

In the derivation of (3.6), we have used the following approximate formulae which are 
valid when <\ x \<^ 1: 

VTTx^l + ^; {l + x)''^l + nx (3.7) 

where in our case x = ~Xexp^a)^" '^^^ derivation of (3.6) shows that (3.5) holds true 
when is much greater than one or when e°- is very large. 

Let [0, to[ be the time interval during which the Lorentzian Hartle-Hawking trajectories 
(hereafter referred to as LHH, see [13]) evolve according to (2.6), (2.7), (3.3) and (3.4). 
We are now interested in deriving solutions of the system 

_dS__dL_ 

Pa — Q — • 

oa oa 

_dS__dL 

in the interval [ta,^f)] during which S can be approximated as in (3.5). Inserting (3.5) into 
this system, and putting N =1, one finds that 

4>= (3-8) 



a 



A/A0t. (3.9) 
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By integrating (3.8) and putting 0^ = (f){ta)i where (f)a = (l>o, we get 



4>{t) 



(4-») 



(2-5)— <'-'•) 



(4-n) 



(3.10) 



for every n 7^ 4. This imphes 



a(t) = a/ exp 



(4-rt) 



{4-n) 



(3.11) 



for every n 7^ 4, where a/ = ^oo^ilita) = 2^exp(ifto). In order to study the case 



(4-n) 



n > 4, we define k = (pa ^ , 6 = (§ - 2)^^, 2; = Therefore (3.11) becomes 



a{t) = a I exp 



b{-z+l) 



{k + bit-ta))^-'+'^ 



(3.12) 



The second term in the curly brackets of (3.12) is 



b{z — 1) n 



(3.13) 



Moreover, at the end of the era during which (3.12) holds, and putting 05 = (f){t},), one has 
from (3.10) that 



th-'ta = 



12 (4-») 



n(n -4)\/A 



(4-») 
2 



]_2 (4-Ti) 



n(n - 4)\/A ' 



(3.14) 



if 0b <^ 0a. Thus the first term in the curly brackets of (3.12) becomes 



(k + hitK - 1 ^^(-^+1) ^ -M 



(3.15) 



From (3.13) and (3.15) we get 



a(tb) = a/ exp 



n 



f-V 



= a/ exp 



30^ 



(3.16) 
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If n = 4, (3.10) and (3.11) are no longer valid; in such a case, the integration of (3.8) 
yields 



(f)= (t)a exp 



-{t-ta) 



(3.17) 



so that (3.9) becomes a — \/A0^ exp(— ^-^(t — ta)), which implies 



a{t) = a I exp 



302 



1 - exp I —{t-ta) 



(3.18) 



Therefore the cosmic scale factor a{ti,) at the end of this era is 



a{ti) = a I exp 



= aj exp 



(3.19) 



if (f)b ^ (f>a, where 



th — ta = 



2VX 



(3.20) 



Furthermore, if n = 3 one finds that 



a{t) = a I exp 



cl>l-{V^a-%t-ta) 



(3.21) 



At the end of the era during which (3.21) holds one has 



tb-ta^4 




1 - 




(3.22) 



Therefore a{tb) = a/exp(0^) if 05 <^ 0a- When ta < t < tb, and t — is very small, 
one has for any n that 



i{t) ^ ai exp{^/X(j)i {t - ta)) = y exp(yA0a' t) 
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which is the correct inflationary formula one would have expected to hold true. It is to be 
emphasized that it would be wrong to put ta = for simplicity in formulae (3.10)-(3.22). 
In fact, t = is the time at which and d vanish according to (2.9), whereas t = ta is 
the time at which S can be approximated by (3.5) up to t — tt,. Therefore during the 
inflationary era and a do not vanish, but can be computed by means of (3.10), (3.11), 
(3.17) and (3.18). This is what also happens in the massive scalar field model [8], in 
which the solutions corresponding to the oscillatory part of the wavefunction start out 
with 0(0) = d(0) = 0, and expand exponentially with = — ^, a = ma | |. 

Finally, when is negative and n is even, (3.5) becomes S ^ . in fact, 

(3.5) is the approximate form of (3.1), which can also be written in the following way: 

1 



3A0" 



e2«A0" - 1 



ye2«A0^ - 1. 



If we choose the square root on the right-hand side to be positive, as we already did in the 
case of positive values of 0, this can only be obtained introducing | | , because otherwise 
0^ could be negative when n is even but ^ is odd. Thus in equation (3.8) is to be 

replaced by the time derivative of | | and in (3.9) 02' is to be replaced by | | ^ . This 
implies that, for any even n, when is negative, | | is a decreasing function of the time 
during the inflationary era. Such a behaviour is typical of a held which starts from 
a large negative value and finally can reach zero before starting to oscillate. Thus the 
relations (3.10)-(3.22) hold true provided that one replaces by | |. 
So far we have derived solutions to the first-order system 

as _ 
da da 

dS__dL 
50 ~ 50 

where 5" is obtained by analytical continuation of the Euclidean action for compact 4- 
metrics and regular matter fields, and L is the Lagrangian defined in (2.1). In particular, 
we have focused our attention on the case when S takes the approximate form (3.5). 
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Some important comments are now in order. In our model the full field equations are the 
second-order differential equations (2.2) and (2.3) plus the Hamiltonian constraint: 

a^^^ -a^ + Xa^(j)'' -1 = 0. (3.23) 

Thus the solution to the full field equations involves three arbitrary constants. On the 
other hand, a wave function of the form Ce*"^ is peaked about the first integral: Pa = ff , 
P0 = This first integral consists of two first-order ordinary differential equations 

(see (3.3) and (3.4)) and so the solution involves just two arbitrary constants. The wave 
function is therefore peaked about a set of solutions which are a subset of the general 
solution. 

Using the Hamilton-Jacobi equation for S, we shall show at first that the solutions of 
the above-mentioned first integral satisfy the full field equations exactly. Finally, we will 
briefly compare the set of solutions picked out by the Hartle-Hawking wave function with 
the general solution. 

The Hamilton-Jacobi equation (2.11) satisfied by S can be cast in the form: 

where g'^^ = diag(— 1, 1), = a, = (f), and V=V(q;, 4>)=\(j)'^ e^"^ — e^'* is the potential 
already defined in section 2. Following Halliwell [21], let us now differentiate (3.24) with 
respect to q^. In so doing we get 

2^-^^ + ^ = 0. (3.25) 
dq^ dqWq^ dq^ 

Let us now introduce, again according to [21], the vector 

^ = 9-^^. (3.26) 
dr oq^ oq^ 



The right choice for the parameter r is its identification with proper time: t = J Ndt. 
In fact one can easily check that in so doing one obtains the well known relation between 
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d, 0, and the momenta Pa and p<^. Therefore ^ = ^ if TV = 1. The insertfon of (3.26) 
into (3.25) yields finally 

and (3.27) is just the form assumed by the field equations. In fact, (2.2) and (2.3) may 
also be cast in the first-order form: 

= -^0"-ie«" (3.28) 

p« = 2e^" - 3A0"e^" (3.29) 
where we have used the following formula for the Hamiltonian: 

H=l{pl-pl)-le''' + ^re'^. (3.30) 

If we now put g'^=0 in (3.27) we obtain (3.28), and if we put =a we obtain (3.29). 

The task remains now to examine the following point: what is so special about the 
solutions picked out by the Hartle-Hawking proposal in comparison to the general solution? 
In fact, if this question is not addressed, it is difficult to see what we have gained by 
computing the wavefunction for our model. 

Indeed, the calculations and the arguments developed so far have shown that the 
solutions about which the Hartle-Hawking wavefunction is peaked are distinguished in 
that they are infiationary. They have this property in view of the fact that the scalar field 
is initially very large so that X(f)'^ acts like an effective cosmological constant in the early 
universe, whereas the initial value of is very small (see (2.9)). However, a member of the 
set of general solutions will not always be infiationary. This is easily understood looking 
more carefully at (2.2) and (3.23). Putting as usual a = ln(a), these relations imply that 

a = -d^ - 202 + xcf)"". (3.31) 



Thus an infiationary solution must be such that d = constant = H>0, d = = 
—H'^ — 20^ + A0". But a member of the set of general solutions might well have an initial 
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value of which is very large. Under such a condition, the right-hand side of (3.31) does 
not vanish if (A^" — d^) <^ 20^. This implies in turn that no inflation is possible if is 
initially so large. 

Therefore the Lorentzian Hartle-Hawking trajectories are very peculiar in that they 
are singularity free in the past (see (2.8)) and inflationary. The full comparison of these 
trajectories with the general solution may be done using the phase-plane method, gener- 
alizing the work done by Belinski et al [22] . This analysis may be an interesting problem 
for further research, but our analysis already shows the main difference between the two 
sets of solutions. 



4. Minimal conditions for a sufficient inflation 



During the inflationary era of the universe, (3.18) becomes 



ait) ^ ai exp[V^0^(t - ta)] = y exp(VA02^). (4.1) 

If we require that the inflationary formula (4.1) satisfies the condition a{t) > aolO^* = 
ao exp (65) in order to solve the horizon and flatness problems [23-25], and if we put t = f3tf), 
where (3 g]0, 1], we find the condition (see (3.20)) 

VA^^,„4^^.og(|)>M!)) (4.2) 

in the case of a A^'* theory. In the same way, we find that the conditions to be satisfied in 

the cases n — 3, n > 4 are, respectively. 



V-X,h.^m(l-j!tl]>J^^±^ (4.3) 



V-,,h^^^m%l>J^^±^. (4.4) 
n(n-4) \(j)bj (3 
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However, a thorough inflationary model also has to solve other problems such as, for 
example, the origin of the energy density fluctuations. Therefore the relations (4.2)- (4.4) 
are just a part of the minimal requirements to be satisfied by our model. 

5. Conclusions 

Many authors [20, 22, 24-28], by using a wide range of techniques, had already considered 
the effect produced in the early universe by a term in the scalar potential of the type 
A0^. In this paper we have studied the general case of a A^" theory in a closed FRW 
minisuperspace model. We have mainly studied the case when (/) is positive and we have 
shown that the case when (f) is negative and n is even can also be taken into account, 
provided that one replaces by | | in all the equations of the theory. But if (p is negative 
and n is odd, this does not give rise to a positive effective cosmological constant in the 
early universe. 

By applying the Hartle-Hawking proposal we have been able to derive the initial 
conditions for the solutions to the Lorentzian field equations (see (2.8) and (2.9)). We 
have also shown that a suitable choice of factor ordering in the Wheeler-DeWitt equation 
(see (2.10)) enables one to extend to our model a result that is already known in the 
case of massive scalar fields: namely, having taken the Euclidean functional integral over 
compact 4-metrics and regular matter fields, the wavefunction is normalizable to one in 
the limit of small 3-geometry and it starts to oscillate when the potential in the Wheeler- 
DeWitt equation is positive and is greater than a number of order one. However, we also 
emphasized that a more careful calculation and a deeper understanding of the semiclassical 
approximation in quantum cosmology are perhaps needed before taking for granted that, 
in general, the wavefunction is a constant on the boundary of minisuperspace. 

Another interesting result is that, with respect to the massive model and for a given 
value of (/), there may be an even larger range of values (i.e. also including smaller values) 
of the cosmic scale factor a{t) for which the potential in the Wheeler-De Witt equation is 
positive. The initial conditions (2.8) are indeed in agreement with this fact: the universe 
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starts off in a non-singular state and the initial value of a{t) is smaller than the one for 
the massive scalar field model, if A is not too small. In the oscillatory region for the 
wavefunction, we have taken the phase S of the WKB approximation to be the analytical 
continuation of the Euclidean action for compact 4-metrics and regular matter fields. The 
LHH trajectories are thus derived by solving the first-order system 



where L is the Lagrangian of the theory (see (2.1)). The solutions of this system are exact 
solutions of the field equations; they have vanishing time derivative in t = and describe 
an exponentially expanding universe at small values of (t — to), where t > ta and [0, ta[ is 
the time interval during which the field (j) has a constant large value. The Hartle-Hawking 
wavefunction is thus peaked about those classical solutions which are inflationary, whereas 
a member of the general set of solutions of the full field equations will not always be 
infiationary. 

The duration of the inflationary era may still be of the order of 10~^^sec or 10~^^sec 
(which is very long if compared to the Planck time) as in the case of the massive scalar 
fleld model, provided that the parameter A is suitably chosen for any value of n. Other 
important conditions to be satisfied are the ones expressed in (4. 2)- (4.4) in order to solve 
the horizon and fiatness problems. Before going over the matter- or radiation-dominated 

phase, the universe may expand up to a maximum value of the order of exp(^^) for any 
n with respect to the value of a{t) at the beginning of the era when exp{—Ht) can be 
neglected in (2.6). As the time passes, decreases, but the approximate formula (3.5) 
from which we have derived most results in section 3 is still approximately valid as far as 
is very large. 

In our opinion it would be particularly interesting to try to solve numerically the field 
equations using the initial conditions derived from the Hartle-Hawking proposal for our 
model, along the lines of the work done by Lafiamme and Shellard [12] for the massive 
model. We think that a thorough examination of this problem, both in the isotropic and 
the anisotropic case, could be of considerable importance in understanding whether or 




dS__dL 
da da 



dS^_dL 
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not the universe described by quantum cosmology will recoUapse to a singularity after a 
maximum expansion. 

Otherwise stated, it may seem quite reasonable that most (or all) minisuperspace toy 
models which are able to drive an inflationary era of the early universe lead to a singularity 
in the future apart from bouncing solutions, but so far this is an unproved conjecture. The 
proof of this result for theories and other models might show an intriguing link between 
inflationary models and the problem of the singularity in the future. In addition, it could 
be even more interesting to show that there are some physically meaningful minisuperspace 
models for which this result does not hold true. 

Finally, we would like to address the attention of the reader to the following point. If 
we multiply (2.1) by —i and if we put N — > —iN, we obtain a Euclidean action Ie which, 
after the substitutions N ^ i, a ^ ia, (f) ^ i(f), becomes 



which is positive definite for all values of n of the form n = 4m, where m = 1, 2, 3, .... This 

seems to suggest that the most interesting X(f)'^ models are just these ones, generalizing 
what was already known in the case of a theory [11, 20]. 
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